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£SJ ■ We present a study of magnetic profiles of umbral dots (UDs) and its consequences 

on the Joule heating mechanisms. Hamedivafa (2003) studied Joule heating using vertical 
component of magnetic field. In this paper UDs magnetic profile has been investigated in- 
cluding the new azimuthal component of magnetic field which might explain the relatively 
larger enhancement of Joule heating causing more brightness near circumference of UD. 



p ' 1. Introduction 

6 

Umbral dots (UDs) are bright features observed in the sunspot umbrae. 
Detailed study of such bright features play a key role in understanding 
the energy transport in sunspots and they exhibit most important physical 
parameters such as temperature, brightness, lifetime, magnetic field, size, 
mass outflows etc. For recent reviews on the subject, see Solanki (2003), 
Thomas and Weiss (2004). 

There are two models for the UDs. First suggested by Parker (1979a, 
b) and then subsequently by Choudhuri (1986). Parker suggested in his 
sunspot model that the region below the visible surface termed as positive 
Wilson depression is made up of individual flaring flux tube embedded in the 
f*"^ \ field free plasma. These tubes merge into an apex like single flux tube just 

above the umbral surface (negative Wilson depression level) of the sunspot. 
The expanding and rising up-flow of plasma creates a gap and develops 
UDs when it reaches to zero Wilson depression surface in sunspot. Further 
Choudhuri (1986) also showed that if the pressure of the plasma plume 
increases at the apex of the static configuration of a field free gas column, 
then it rises to certain height where the magnetic field pressure suddenly 
becomes negligible. Finally the trapped field free gas bursts in to a column at 
a speed of about 10 km s~ 1 forming a UD. In another magnetohydrodynamic 
model, UDs are considered as tops of convective cells present in homogeneous 
magnetic field (Knobloch and Weiss, 1984). Later Degenhardt and Lites 
(1993) speculated that UDs are thin vertical magnetic tubes with a reduced 
magnetic field strength, a temperature enhancement and material upward 
flow from the bottom, embedded in a stationary sunspot umbra. 
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Some experimental studies have revealed that the lifetime of UDs is in 
the range from 15 min to 2 hrs with average diameter of 150-300 km and 
the relative intensities with regard to umbra of sunspot vary from 1.1-2.6 
(Lites et al, 1991; Ewell, 1992; Sobotka et al, 1997a; Tritschler and Schmidt, 
2002; Sobotka and Hanslmeier, 2005). The number of investigators have also 
observed mass upward motion within the range of 0.3-1 km s _1 (Pahlke and 
Wiehr, 1990; Lites et al, 1991; Wiehr, 1994; Rimmele, 1997, 2004). 

The role of magnetic profile has been unique for some specific properties 
of UDs and in one of the studies of magnetic profiles, Choudhuri (1986) 
established that there has been a reduced magnetic strength in the UD col- 
umn relative to the surrounding umbra. Some spectroscopic observations of 
Adjabshirzadeh and Koutchmy (1983), Pahlke and Wiehr (1990), and Wiehr 
and Degenhardt (1993) have shown that positive magnetic field gradients in 
umbral dots are of about 20%, whereas in the high resolution study of Lites et 
al. (1991) have revealed that magnetic strength is not significantly different 
from that the surrounding darker portions of the umbrae. Also recent study 
(Socas-Navarro et al., 2004) showed that there is more inclined fields in 
UDs and the vertical gradient of field may have opposite signs in the UD 
and dark background. Since all observations are obtained often with a low 
spatial resolution restricted for a small area so that the complete magnetic 
vector remains unknown. Thus these observational results point out that a 
good knowledge of magnetic vector field is crucially required. However, in 
recent simulation study of Schiissler and Vogler (2006), they have manifested 
that nearly field free upflow plumes and the UDs are a natural result of 
convection in a strong initially monolithic magnetic field. 

In the subsequent Section 2 magnetic field profile for UDs has been 
discussed whereas in Subsection 2.2, we have tried to develop the modified 
magnetic field profile by introducing a proposed azimuthal component B^r) 
crucially required with vertical component B z (r) and its effect on Joule 
heating. The resultant magnetic vector obtained for the UDs might be at- 
tributed to the increased current density which probably may enhance Joule 
heating power causing a relatively more brightness near the circumference 
as compared to the center of the UDs. 



2. Magnetic Profile of Umbral Dots and Joule Heating 

We still know very little about the structure and the nature of the magnetic 
field of UD, as there are few high-resolution evidences available (Lites et al., 
1991; Socas-Navarro et al, 2004) for the vector magnetic field that can pro- 
vide correct magnetic field structure in UDs. In spite of these observational 
challenges, the right choice of magnetic profile may produce better consis- 
tency with observations. In the next Subsection 2.1 we first give an overview 
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of vertical magnetic field and then we propose an additional azimuthal 
magnetic field with the vertical component in the subsequent Subsection 
2.2. 

2.1. Vertical component of magnetic field in umbral dots 

As a matter of fact Joule heating power is partially responsible for brightness 
of the UDs and is directly governed by current density, which in principle is 
attributed to magnetic field profile (Garcia de la Rosa, 1987). Thus it was 
assumed for simplicity that magnetic field vector has only vertical paral- 
lel component in the UD column (Hamedivafa, 2003) and is a function of 
distance from axis of the UD column as 

B = B z (r)8L z (1) 

where a z is the unit vector along z-axis normal to the photosphere of the Sun 
and r is the radial component in cylindrical coordinates (r, (ft, z). Hamedivafa 
(2003) investigated Joule heating as brightening mechanism for umbral dots 
and he assumed one of the magnetic profiles given as 

B z (r) = B - 1 B exp^- (2) 

where Bq is saturated magnetic field in umbra of the sunspot, 7 is fraction 
of the field redution and u is the maximum radius of the UD column at any 
instant. Hamedivafa and Sobotka (2004) found direct observational evidence 
for Joule heating in some of the UDs. 

2.2. Azimuthal component of magnetic field in UDs 

A prudential overview of the existing models for the magnetic field of UDs 
is required since the magnetic nature of bright features in sunspot umbrae 
is not fully understood yet. The observational results of Lites et al. (1991) 
pointed out that UDs magnetic field strength is not very different from 
umbral field and also Socas-Navarro et al. (2004) showed that UDs have 
inclined fields. This promulgated the following two major assumptions which 
have been contemplated in theoretical investigations: 

1. B z (r) is the force free axial vector field for the expanding and rising 
parallel axial flow of plasma plume in UD column. It develops cusp- 
like shape due to steeper pressure gradient of piled-up plasma in UDs 
leading to a drastic decrease of magnetic field strength in the upper layer 
of plumes than the proximate circumference of UD. 

2. Steady flow of material in UD column does not give a temporally chang- 
ing electric field i.e. there is no changing electric flux. 
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These assumptions deduce a basic magneto-hydrodynamic equation 

(V x B) 2 oc fi J z (3) 

which is valid for azimuthal component. Where is the magnetic perme- 
ability across the active region and the current density J z is associated with 
a convectively unstable and oscillating vertical slab of plasma, sandwiched 
between the regions of B z (r) with vertical wavelength of 100 km and of 
period 100 s (Parker, 1979b). Moreover it may be expressed as a possible 
electric current parallel to the force free B z {r) within r < u column as 

/ = j> J z -da (4) 

Here r is the radial distance from the UD axis; u is the maximum radius of 
UD at which an undisturbed magnetic field strength of umbra exists. The 
high temperature of the hot unstable outflowing material would be sufficient 
to yield an ionized form, creating enough electrical current in these column. 
These concepts led us to believe that there should be azimuthal magnetic 
component B,p(r) attached to the thin UD column, due to electric current 
of material at high temperature. Therefore, the magneto-hydrodynamic cal- 
culations under some boundary conditions for evolution of these flux tubes, 
diagnose UDs of enhanced temperature and high intensity relative to sur- 
rounding umbra. However this current would be able to generate azimuthal 
magnetic field B^r) with the help of Ampere circuital law, Equation (3), 
and B(f,{r) may be expressed as 

7/i/r 

B(f>(r) = - — j when r < u (5) 
7 fj/I 

BAr) = when r > u. (6) 

9y ' 2irr w 

Here 7 is the fraction of the field strength reduction on the central axis of 

the UD column with range 1 > 7 > 0. Hamedivafa (2003) has revealed that 

the magnetic field reduction is not beyond the radius of the bright UDs but 

saturates after r/u ~ 2. The contribution of B^r) is linear from the center 

to the circumference of the UD and then hyperbolically decreases beyond the 

UD column. The inclusion of B^r) in the proposed model may be justified 

on the basis of following assumptions: 

1. The center of the azimuthal component coincides with the symmetric 
axis of the UDs. It has been assumed to resolve the ambiguity between 
center of B^{r) and vertical component B z (r) which yields more or 
less fine radial magnetic structure inside UDs producing a temperature 
stratification and a height independent values for magnetic field strength 
with respect to line of sight. 
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Figure 1. Normalized magnetic field strength of UD versus normalized distance from the 
UD centre. 

2. Magnetic field of the UDs consists of the two components, -B^(r) and 
B z {r) to provide a consistent feature of magnetic profile in such a way 
that the magnetic field inside an UD is weaker than surrounding field 
and large magnetic field gradient is present at the boundary of UD. The 
component B^r) might contribute more power to Joule heating. 

3. B^ir) at the circunference of an UD becomes comparable to the observ- 
able magnetic field, Bq of umbrae. 

Thus resultant B(r) may be written as 

B(r) = B*(r)a* + £*(r)a*. (7) 

Normalised B(r), B^r) and B z (r) with Bq are plotted as shown in Figure 
1. The effective current density in UDs is calculated by the magnetic field 
as 

Jt = z Vx b M ( 8 ) 
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or the value of Jt is calculated in cylindrical co-ordinate system by the 
followig determinant 



A _d_ A 

dr rd(j> dz 

S 5 - 7 5oe X p^i 



— c — r 2 \ 2r 

— 7-B exp x — 



CJ/J.I 

8tt 2 u 2 



Now the Joule heating power can be given as: 

4tt /-oo 



Q' = — v4(^r)dr, 
c z Jo 



where r\ is electrical resistivity. 



Jt • Jt 



4r 2 



16vr 2 



-2r 2 



+ 



c 2 7 V/ 2 
64vr 4 u 4 



(9) 



(10) 



(11) 



From Equation (10) and (11), we get 



Q' 



2 E>2 



2t ?7 2 j B ( 



u 4 \ rij 2 /j, 2 I 2 u 2 
Y) + 4(4tt 2 u 2 ) 



u' 1 



(12) 



We assume that when r — > u then \B,p(r)\ — ► |i?o(^)| , i.e. the maximum 
strength of azimuthal component at the boundary of UD is comparable 
with the undisturbed magnetic field Bq of the umbra outside of the UD. 
Thus from Equation (5) we get 



Bn 



^filu 7 //I 



or 



Let 



2ttu 2 2iru 
2 (2— - x2 

~2 



2 2 _ B 2 (2nuy 

7' 



(13) 
(14) 



U 1 b 



8vr ' 



(15) 



With the help of Equations (12) - (15), the final power can be given as 



Q' 



*+4 

7 . 



(16) 



The Joule Heating power is proportional to 2irr]UyB whereas Q = 2ttt]U 1 b 
is calculated by Hamedivafa (2003). The factor in the square bracket rep- 
resents the enhancement in the power due to the current generated from 
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the magnetic field at r close to u. For a smaller value of 7, larger would be 
the Joule heating power. This justifies the special features observed in the 
experimental studies of Lites et al. (1991) that the field strength within the 
umbrae vary at large scale (1400 - 2400 Gauss) and this large scale variation 
of field being inversely correlated with the umbral dot intensity. 

3. Discussion and Conclusions 

The result presented in this report depends upon the following crucial as- 
sumptions: 

1. The B^ir) is developed just at the middle of UDs and increases linearly 
with r inside the UD diameter up to the maximum value B . It decreases 
hyperbolically outside. 

2. The current density Jt is uniform within the UDs. 

3. Magnetic profile of UDs consists of at least two magnetic components, 
B^{r) and B z (r) which stem not necessarily from the same geometrical 
height. 

We conclude enhanced field on peripheral surface of UD as shown in 
Figure 1. The large peripheral flux is attributed to the additional azimuthal 
magnetic field component, which is produced by axial current / due to in- 
trusion of hot plasma. Thus modified UD model obtained by adding B^ij) 
into B z {r) helps to explain variation in brightness of the UDs. The relative 
Joule heating 

AQ = f = l + l (17) 

Table I. The relative Joule heating with diffrent values of 7 

7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
A Q 101.0 26.0 12.11 7.25 5.00 3.78 3.04 2.56 2.23 2.00 



Fractional Joule heating is much more at lower values of 7 but at higher 
values decreases slowly. The estimated Joule heating power is expected 
to justify the inverse correlation of brightness and magnetic field gradient 
(Lites et at., 1991) and the resultant magnetic vector supports more inclined 
magnetic field vector (Socas-Navarro et al., 2004) of the UD. 

Moreover, these conclusions can be justified by high resolution spec- 
tropolarimetric data from ground based instrument such as Diffraction Lim- 
ited Spectropolarimeter (DLSP) (Sankarsubramanian et al, 2004) at Dunn 
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Solar telescope (NSO) and the spectropolarimeter onboard the satellite 
Hinode. 
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